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Abstract 

The first observation of the decay B° — > D°K + K~ is reported from an analysis of 
0.62 fb _1 of pp collision data collected with the LHCb detector. Its branching frac- 
tion is measured relative to that of the topologically similar decay B° — > D°it + it~ 

to be _ 

B(B° -»• D°K+K-) 

\ =- ^ = 0.056 ± 0.011 ± 0.007 , 

B (B° -»■ D°7r+7r-) 

where the first uncertainty is statistical and the second is systematic. The signifi- 
cance of the signal is 5.8 a. Evidence, with 3.8 a significance, for B® — > D°K + K~ 
decays is also presented. The relative branching fraction is measured to be 

B (B° s D°K + K- 



B (B° -»■ D°K+K~ 



0.90 ±0.27 ±0.20. 



These channels are of interest to study the mechanisms behind hadronic B decays, 
and open new possibilities for CP violation analyses with larger data sets. 
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The precise measurement of the angle 7 of the CKM Unitarity Triangle [1]|2] is one 
of the primary objectives of flavour physics experiments. Prior to the start of LHC data- 
taking, the combination of measurements with the decay mode B + — > DK + , where D 
denotes a neutral charmed meson that is an admixture of D° and D°, gave a constraint 
on 7 with an uncertainty of around 20° [3j. Recent results from LHCb on B + — > DK + |I] 
have helped to reduce this uncertainty, but the use of additional channels to improve 
further the precision is of great interest. The as- yet unobserved decay B® — > Dip is one 
of the modes with potential to make a significant impact on the overall determination 
of 7 j5]-[7j. Moreover, a Dalitz plot analysis of B° s — > DK + K~ can further improve the 
sensitivity to 7 due to heightened sensitivity to interference effects, as well as allowing a 
determination of 4> s , the CP- violating phase in the B®-B® system, with minimal theoret- 
ical uncertainties |8]. 

The first step in the programme towards the measurement of 7 using the B® — > 
DK + K~ decay is the observation of the channel. In this Letter the results of a search 
for neutral B meson decays to D°K + K~ are presented. The quantities measured include 
small contributions from decays to D°K + K~ . The inclusion of charge conjugate modes 
is implied throughout. 

The analysis uses 0.62 fb _1 of LHC collision data at a centre-of-mass energy of 7TeV 
collected with the LHCb detector during 2011. In high energy pp collisions all b hadron 
species are produced, so both B° and B® decays are searched for simultaneously. The 



ll These 
in this 



decay B° — > D°K + K~ can be mediated by the decay diagrams shown in Fig. 
are a ly-exchange diagram similar to that for the decay B° — > D~K + (9|[To 
case an excited state that decays to D°K~, such as D*^"(2573), would be produced), 
and a colour-suppressed tree diagram producing D°h°, where h° is a light unflavoured 
meson such as a (980) that subsequently decays to K + K~. Rela ted B decays with ss 
production, B+ D°K + K^° [11] and B + Di* y K + n 
to have branching fractions of 0(1O -4 ). 
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have been measured 




Figure 1: Sample decay diagrams that contribute to the B° — > D°K + K final state via 
(left) ^-exchange, (right) rescattering from a colour-suppressed decay. 



The LHCb detector 14 is a single-arm forward spectrometer covering the 
pseudorapidity range 2 < r\ < 5, designed for the study of particles containing b or c 
quarks. The detector includes a high precision tracking system consisting of a silicon-strip 
vertex detector surrounding the pp interaction region, a large-area silicon-strip detector 
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located upstream of a dipole magnet with a bending power of about 4 Tm, and three sta- 
tions of silicon-strip detectors and straw drift tubes placed downstream. The combined 
tracking system has momentum resolution Ap/p that varies from 0.4% at 5 GeV/c to 0.6% 
at 100 GeV/c, and impact parameter (IP) resolution of 20 |xm for tracks with high trans- 
verse momentum (pr)- Charged hadrons are identified using two ring-imaging Cherenkov 
(RICH) detectors. Photon, electron and hadron candidates are identified by a calorime- 
ter system consisting of scintillating-pad and pre-shower detectors, an electromagnetic 
calorimeter and a hadronic calorimeter. Muons are identified by a system composed of 
alternating layers of iron and multiwire proportional chambers. The trigger consists of a 
hardware stage, based on information from the calorimeter and muon systems, followed 
by a software stage which applies a full event reconstruction. In this analysis, signal can- 
didates are accepted if one of the final state particles created a cluster in the calorimeter 
with sufficient transverse energy to fire the hardware trigger. Events that are triggered 
at hardware level by the decay products of the other b hadron in the pp — > bbX event are 
also retained. 

The software trigger requires characteristic signatures of b hadron decays: at least 
one track, with high p^ and a large IP with respect to any primary interaction vertex 
(PV) [l5], that subsequently forms part of a two-, three- or four-track secondary vertex 
with a high sum of the pi of the tracks and significant displacement from the PV [16) . 
In the offline analysis, the software trigger decision is required to be due to the candidate 
signal decay. 

Candidates that are consistent with the decay chain B9s — > D°K + K~ , D° — > K + n~ 
are selected. In order to reduce systematic uncertainties in the measurement, the topo- 



logically similar decay D°7t + tt~, which has previously been well studied 17,18 , is used 
as a normalisation channel. The D° candidate invariant mass is required to satisfy 
1844 < rriKn < 1884MeV/c 2 . Tracks are required to be consistent with either the kaon or 
pion hypothesis, as appropriate, based on particle identification (PID) information from 
the RICH detectors. All other selection criteria were tuned on the D°7t + tt~ channel. The 
large yields available in the normalisation sample allow the selection to be based on data, 
though the efficiencies are determined using Monte Carlo (MC) simulated events. For the 
simulation, pp collisions are generated using Pythia 6.4 [19] with a specific LHCb config- 
uration 20 . Decays of hadronic particles are described by EvtGen 21 . The interaction 



of the generated particles with the detector and its response are implemented using the 
Geant4 toolkit [22) as described in Ref. [23 . 



Selection requirements are applied to obtain a clear signal peak in the D°7t + tt~ normal- 
isation channel. The selection includes criteria on the track quality of the tracks forming 
the signal candidate, their p, p^ and inconsistency with originating from the PV (xip)- 
Requirements are also placed on the corresponding variables for candidate composite par- 
ticles (D°, B9s) together with restrictions on the consistency of the decay fit (Xvertex)> the 
flight distance significance (xlight)' an d the angle between the momentum vector and the 
line joining the PV to the B? s) vertex (cos 6^) [24]. 



Further discrimination between signal and background categories is achieved by cal- 



culating weights for the remaining D n + n candidates 25 . The weights are used by the 
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NeuroBayes neural network package [26] to maximise the separation between categories. 
A total of 15 variables are used in the network. They include the Xip °f the four candi- 
date tracks, the Xipi X V ertex> Xfiight an d cos^dir of the D° and B^ candidates, and the B9 g \ 
candidate pp- Variables describing the pt asymmetry and track multiplicity in a 1.5 rad 
cone j4j around the candidate flight direction are also used. The input quantities to 
the neural network only depend weakly on the kinematics of the B^ decay. A require- 
ment on the network output is imposed that reduces the combinatorial background by an 
order of magnitude while retaining about 80 % of the signal. No bias is observed by using 
data driven selection requirements. 

To improve the B9\ candidate invariant mass resolution, the four-momenta of the 

tracks from the D° candidate are adjusted so that their combined invariant mass matches 
the world average value (3). An additional B° mass constraint is applied in the calculation 
of the Dalitz plot coordinates, which are used in the determination of event-by-event 
efficiencies. A small fraction (~ 5 % within the mass range described below) of candidates 
with invariant masses far from the B^ peak fail the mass constrained fit, and are removed 
from the analysis. 

To remove a large potential background from B° — > D*~ (2010)7r + , candidates in the 
D°7r + 7r~ sample are rejected if m^-m^ (for either pion charge) lies within ±2.5MeV/c 2 
of the nominal D*~-D° mass difference (3). Candidates in the D°K + K~ sample are 
also rejected if the invariant mass difference calculated under the pion mass hypothesis 
satisfies the same criterion. This removes 3.3% of D°K + K~ candidates. Less than 1% 
of D°K + K~ combinations are rejected by requiring that the pion from the D° candidate 
together with the two kaons do not form an invariant mass in the range 1950-1975 MeV/c 2 , 
which removes potential background from B® — > DJK ± decays. 

After all selection requirements are applied, less than 1 % of events with at least one 
candidate also contain a second candidate. Such multiple candidates are retained and 
treated the same as other candidates; the associated systematic uncertainty is negligible. 

In addition to combinatorial background, candidates may be formed from misidenti- 
fied or partially reconstructed B9-. decays, or from B ^ decays to identical final states 
but without intermediate charmed mesons (referred to below as charmless peaking back- 
ground). Contributions from partially reconstructed decays are reduced by requiring the 
invariant mass of the B9s candidate to be above 5150 MeV/c 2 . Sources of misidenti- 
fied backgrounds are investigated using simulation. Most potential sources are found to 
have a broad invariant mass distribution, and are absorbed in the combinatorial back- 
ground shape used in the fit described below. Backgrounds from A b — > D°pK + and 
A° b — > D°pn + E7J, B° — > D°K + 7i~ and B® — > D°K~n + decays may, however, give contri- 
butions with distinctive shapes and therefore need to be included in the fit. 

The contributions from charmless peaking background are investigated using candi- 
dates, reconstructed without the D° mass constraint, in sideband regions around the 
D° mass. The distributions are fitted with double Gaussian signal and linear back- 
ground probability density functions (PDFs). Extrapolating to the D mass signal region, 
773 ± 30 (126 ± 18) charmless background decays are expected in the B° — > D°7t + tt~ 
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(B° — > D°K + K~) distributions. No peaking background is observed in the B° s region. 

The signal yields are obtained from unbinned maximum likelihood fits to the D 7t + tt~ 
and D°K + K~ invariant mass distributions in the range 5150-5600 MeV/c 2 . There are 
14214 D°ir + ir- and 2990 D°K + K~ candidates. The D^+n' fit includes a double Gaus- 
sian shape for signal, together with an exponential component for partially reconstructed 
background, and a PDF for A b — > D°X decays modelled using a non-parametric func- 
tion obtained from simulation. The D°K + K~ fit includes a second double Gaussian 
component to account for the possible presence of both B° and B® decays, and peaking 
background PDFs for A° b -»■ D°pK + , B° -»■ D°K+ir- and B° s -»■ D K-tt + , all modelled 
using non-parametric functions. The shape of the combinatorial background is essentially 
linear, but is multiplied by a function that accounts for the fact that candidates with high 
invariant masses are more likely to fail the B?^ mass constrained fit. 

The result of the fit to D°n + ix~ candidates is shown in Fig. [2J There are nine free 
parameters in this fit: the double Gaussian peak position, core width and fraction in 
the core, the linear slope of the combinatorial background and the exponential shape 
parameter of the partially reconstructed background, and the yields of the four categories. 
The relative width of the broader to the core Gaussian component is constrained within 
uncertainty to the value obtained in simulation. The fit yields 8060 ± 150 B° — > D°7t + 7t~ 
decays, including charmless peaking background. 




Figure 2: Fits to the B9^ candidate invariant mass distributions for the (left) D°-k + 7t 
and (right) D°K + K~ samples. Data points are shown in black, the full fitted PDFs as 
solid blue lines and the components as detailed in the legends. Yields of the partially 
reconstructed and peaking backgrounds are all small for the D°K + K~ sample. 



Since the fit to D°K + K~ candidates, shown in Fig. [2], has more components, additional 
constraints are imposed in order to improve the stability of the results. The parameters 
of the double Gaussian shapes are constrained to be identical for B° and B° s signals, with 
an offset in their mean values fixed to the known B°-B° s mass difference j3| . The slope of 
the combinatorial component is constrained to the value obtained in the fit to D° mass 
sideband events. The exponential shape parameter is constrained to the value obtained 
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in the D 7i + 7T^ fit. The fit yields 558 ± 49 B° — > D°K + K~ decays, including charmless 
peaking background, and 104 ± 29 — > D°K + K~ decays. All background yields are 
consistent with their expectations within uncertainties. 

The ratio of branching fractions is obtained after subtracting the charmless peak- 
ing background, and applying event-by-event efficiencies as a function of the Dalitz plot 
position 



B (5° D°K + K~) Q ^ _ N cmr (D°K + K~) (l 



NP e ^(D°K+K-) 



TO») = ^^ -v r^-'/ , (1) 



v ' y iV (_L> u 7r+7r ) 

where N is the yield obtained from the fit, iV peak is the charmless peaking background 
contribution, and the efficiency corrected yield jV corr = £\ Wi/e\ ot . Here the index i runs 



over all candidates in the fit range, Wi is the signal weight for candidate i 25 from the fit 
shown in Fig. 2 and e\ ot is the efficiency for candidate i, which depends only on its Dalitz 
plot position. The statistical uncertainty on the branching fraction ratio incorporates 
the effects of the shape parameters that are allowed to vary in the fit, the dilution due 
to event weighting, and the charmless peaking background subtraction. Most potential 
systematic effects cancel in the ratio. 

The PID efficiency is measured using a control sample of D*~ — > D°7T~ , D° — > K + n~ 
decays to obtain background-subtracted efficiency tables for kaons and pions as functions 
of their p and pr [28] . The kinematic properties of the tracks in signal decays are obtained 
from simulation, allowing the PID efficiency for each event to be obtained from the tables 
taking into account the correlation between the p and px values of the two tracks. The 
other contributions to the efficiency (detector acceptance, selection criteria and trigger 
effects) are determined from simulation, and validated using data. All are found to be 
approximately constant across the Dalitz plane, apart from some modulations seen near 
the kinematic boundaries. 

The Dalitz plot distributions obtained from the signal weights are shown in Fig. |3j 
The B° — > D°7i + tt~ distribution shows contributions from the p°(770) and /2(1270) reso- 
nances (upper diagonal edge of the Dalitz plot) and from the DT~ (2460) state (horizontal 



band), as expected from previous studies of this decay 17, 18 . The B° — > D°K + K~ 
distribution shows a possible contribution from the _D*2~(2573) resonance, together with 
an enhancement of events at low K + K~ invariant mass (upper diagonal edge). 

The branching fraction of the B® decay to D°K + K~ is measured relative to that of 
B° to the same final state. Due to the low yield in this decay, an event-by-event efficiency 
correction is not used. The ratio of branching fractions is instead determined as 

B (B° D°K + K~) _ R{B o B o, _ ffsY 1 N(B° -> DKK) 



B(B° ^ D°K+K-) v s ' ' \f d J N(B° -> DKK) - N^(B° ->■ DKK)' 

_ (2) 

The ratio of fragmentation fractions is f s / fd — 0.267 ioirao 
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Systematic uncertainties are assigned to both branching fraction ratios due to the 
following sources. The variation of efficiency across the Dalitz plot may not be correctly 
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Figure 3: Dalitz plot distributions for (left) B° ->■ D ^^ and (right) B° ->■ D°K + K~ 
obtained from the signal weights. Note that these distributions contain contributions 
from charmless peaking backgrounds. 



modelled in simulation. The difference, 6.7%, between the nominal result for R(B°,B°) 
and that obtained using Dalitz plot averaged efficiencies is conservatively taken as an 
estimate of the associated systematic uncertainty. The fit model is varied by scaling 
the B®/B° PDF width ratio to account for their different masses, removing components 
with small yields, adding components for potential background from B° s — > D*°K*° and 
B° s — > D*°K + K~ , and varying the linear parameter of the combinatorial background PDF 
within uncertainties from the fit to the D° sidebands used to estimate the charmless peak- 
ing background. Together these contribute 10.7% (19.9 %) to R(B°, B°) (R(B° S , B )). An 
uncertainty of 1.5% is assigned due to the charmless peaking background subtraction 
procedure. Possible biases in the determination of the fit parameters are investigated 
using MC pseudoexperiments, leading to 1.5% (3.4%) uncertainty on the R(B°,B°) 
(R(BlB°)). 

In addition, the possible differences in the data/MC ratios of trigger and PID efficien- 
cies between the two channels (both 2.0%) and the effect of the D+ veto (1.7%) affect 
only R(B°,B°). The uncertainty on the quantity f s /f d (7.9%) affects only R(B° S ,B°). 
The total systematic uncertainties are obtained as the quadratic sums of all contributions. 

A number of cross-checks are performed to test the stability of the result. The data 
sample is divided by dipole magnet polarity, data taking period and trigger category. 
Candidates were divided based upon the hardware trigger decision into three groups; 
events in which a particle from the signal decay created a large enough cluster in the 
calorimeter to fire the trigger, events that were triggered independently of the signal 
decay and those events that were triggered by both the signal decay and the rest of the 
event. The neural network and PID requirements are tightened and loosened. The PID 
efficiency is evaluated using the kinematic proprties from D n + n^ data instead of from 
simulation. The charmless peaking background contribution is determined from the upper 
and lower D° mass sidebands separately. All give consistent results. 
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The significances of the signals are obtained from the changes in likelihood in fits to 
data with and without signal components, after accounting for systematic uncertainties 
and for charmless peaking background in B° — > D°K + K~ only. They are found to be 
5.8(7 and 3.8 a for B° ->- D°K + K-_and B° s ->- D°K + K~ respectively. 

In summary, the decay B° — > D°K + K~ has been observed for the first time, and its 
branching fraction relative to that of B° — > D°7t + 7t~ is measured to be 

B(B° ^ D°K+K-) 

— ^ = ^ = 0.056 ± 0.011 ± 0.007 , 

B(B° D°n+n-) 

where the first uncertainty is statistical and the second is systematic. Using the known 
value of B (B° DV^r") = (8.4 ± 0.4 ± 0.8) x 10" 4 [Tt] , this gives 

B (B° D°K + K~) = (4.7 ± 0.9 ± 0.6 ± 0.5) x 10~ 5 , 



where the third uncertainty arises from B — > D°n + n ) . Evidence for the B 
D°K + K~ decay has also been found, with relative branching fraction 

B (B° s D°K + K~) 



o 



B (B° D°K+K~) 



0.90 ±0.27 ±0.20. 



A future study of the Dalitz plot distributions of these decays will provide insight into 
the dynamics of hadronic B decays. In addition, the B® — > D°K + K~ decay may be used 
to measure the CP violating phase 7. 
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